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The vertical structure of Jupiter's equatorial zonal wind above 
the cloud deck, derived using mesoscale gravity waves 

C. Watkins 1 and J. Y-K. Cho 1 
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Data from the Galileo Probe, collected during its de- 
scent into Jupiter's atmosphere, is used to obtain a ver- 
tical profile of the zonal wind from ~0.5 bar (upper tro- 
posphere) to ~0.1/ibar (lower thermosphere) at the probe 
entry site. This is accomplished by constructing a map of 
gravity wave Lomb-Scargle periodograms as a function of 
altitude. The profile obtained from the map indicates that 
the wind speed above the visible cloud deck increases with 
height to ~150 ms _1 and then levels off at this value over 
a broad altitude range. The location of the turbopause, 
as a region of wide wave spectrum, is also identified from 
the map. In addition, a cross-equatorial oscillation of a jet, 
which has previously been linked to the quasi-quadrennial 
oscillation in the stratosphere, is suggested by the profile. 



1. Introduction 

The horizontal structure of Jupiter's zonal (east-west) 
winds at the level of the visible cloud deck on Jupiter have 
been observed for many decades [e.g., Limaye, 1986; In- 
gersoll, 2004; Garcia-Melendo and Sanchez Lavega, 2001; 
Vasavada and Showman, 2005; Li et ai, 2006]. In contrast, 
the vertical structure of zonal winds is not well known - 
especially away from the cloud deck and in the equatorial 
region. This is because of the lack of usable tracers and 
the breakdown of the thermal wind relation near the equa- 
tor, preventing latitudinal temperature measurements to be 
related to the vertical wind shear. 

In December 1995 the Galileo probe entered Jupiter's up- 
per atmosphere at 6° N, above a 5/im hot spot in the north 
equatorial belt. During the descent the Atmospheric Struc- 
ture Instrument on board collected information about the 
density, pressure, and temperature of the atmosphere [Seiff 
et at, 1998]. Analysis of the temperature profile in the ther- 
mospheric and stratospheric regions identified discernible 
perturbations which have been interpreted as manifestations 
of internal, or vertically propagating, gravity waves [Young 
et al, 1997, 2005; Matcheva and Strobel, 1999]. Vertical tem- 
perature profiles for Jupiter's atmosphere have also been ob- 
tained from occultation studies using stars and other space- 
craft. A number of these profiles contain oscillations that 
have been characterised as manifestations of internal grav- 
ity waves as well [French and Gierasch, 1974; Lindal, 1992; 
Hubbard et al, 1995; Raynaud et al, 2003, 2004]. 

Gravity waves are oscillations of fluid parcels about their 
altitudes of neutral buoyancy [Gossard, 1975]. The waves 
are a common feature of stably stratified atmospheres. They 
have been captured in images of Jupiter's clouds, typically 
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near the planet's equator [Flasar and Gierasch, 1986; Ar- 
regi et al, 2009; Reuter et al, 2007]. Fig. 1 shows an 
example. Internal gravity waves (hereafter simply gravity 
waves) grow in amplitude, due to the fall in background at- 
mospheric density. The dynamics of such waves is described 
by the Taylor-Goldstein equation [Taylor, 1931; Goldstein, 
1931], 



d 2 w . 
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Here w is the perturbation in the vertical wind, which is 
adjusted for the amplitude growth and also assumed to be 
oscillatory (wave-like) in the horizontal direction and time; 
F is a function that represents the source and dissipation of 
waves; and, m is the vertical wavenumber, which depends 
on the physical properties of the medium in which the waves 
propagate. Specifically, 
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where N is the Brunt -Vaisala (buoyancy) frequency; / is 
the intrinsic phase speed, c — uo, where c is the horizon- 
tal phase speed and Mo is the zonal wind; H is the density 
scale height; and, k is the horizontal wavenumber. In gen- 
eral, all of the variables depend on altitude z. Crucially, 
the zonal wind profile wo(^) can be obtained by solving for 
I = I(H, k, N, m) in equation (2). In this letter we present 
a profile from this inversion that span much larger range of 
altitude above the 1 bar level (range of ~ 500 km) than in 
past analyses. 

2. Method 

For the inversion we compute H from the collected den- 
sity profile. We also use k = 2tt/300 km" 1 fr om the average 
value obtained in analysis of Voyager images [Flasar and 
Gierasch, 1986]. Other analyses have reported slightly dif- 
ferent values (e.g., average of 2-7r/165 km -1 by Arregi et al. 
[2009]); however, as long as k 2 <C m 2 (as is the case here), 
the value of k does not significantly alter the result of the 
inversion. N and m are obtained from the potential tem- 
perature, 6 = T(p le {/p) K , related to specific entropy of the 
atmosphere; here T is the temperature, p is the pressure, 
p rc f = 1 bar is a constant reference pressure and k = R/c p 
with R the specific gas constant and c p the specific heat at 
constant pressure. To obtain N and m, we decompose 6 into 
a mean background 6 (extracted using a spatially-moving fit- 
ting window) and a small perturbation A8 about the mean 
(presumed caused by mesoscale gravity waves). Such decom- 
position of temperature is standard in gravity wave studies 
[Lindzen, 1990; Nappo, 2002; Young et al, 1997], and is rea- 
sonable here given that the spatial scales of the background 
and waves are well separated. Further, the horizontal dis- 
tance travelled by the probe in the region studied, 3300 km, 
is small compared to Jupiter's circumference and is com- 
pletely within the hot-spot entered by probe [Orton et al, 



1 



X - 2 



WATKINS AND CHO: JUPITER'S EQUATORIAL ZONAL WIND 



1998]; hence, zonal variation in the background over this 
distance is not expected to be significant. 

Fig. 2a shows the relative perturbation, A9/9, result- 
ing from a 75 km- wide moving window. We have checked 
that the result is not affected by the choice of window size, 
by varying the size from 55 km to 85 km. Fig. 2b shows 
N = [g d(ln ^/dz] 1 / 2 resulting from the decomposition; here 
g — g(z) is the gravity. In the figure note that N is much 
smaller in the region below ~20 km than that above (the 
stratosphere) — in agreement with Magalhaes et al. [2002], 
whose observations give N « 6 x 10 s -1 in the lower re- 
gion. This represents a possible ducting region, a source for 
the waves in the probe data. 

Although w information was not collected by the probe, 
m(z) can still be obtained from A8 through the polarization 
relation [e.g., Watkins and Cho, 2010]. For this we gener- 
ate a series of Lomb-Scargle periodograms [Scargle, 1982] 
of the A8 data, which is non-uniformly spaced in z. The 
non-uniform spacing renders analysis by standard Fourier 
or wavelet transforms unsuitable. One periodogram at each 
z is generated using a smoothing window. We have also 
verified in this procedure that the obtained result varies 
little between different sized smoothing windows. All the 
periodograms are subsequently combined to produce a two- 
dimensional map of the wave spectral energy density as a 
function of wavenumber and altitude, £ — £(m*,z), shown 
in Fig. 3; here m* is the vertical wavenumber m prior to 
an adjustment for wave propagation geometry. Before dis- 
cussing this adjustment, we discuss three features which are 
already apparent in the £ map. 

3. Results 

First, gravity waves identified previously [Young et al., 
2005] are recovered (Wl and W2 in the map). This agree- 
ment gives confidence in our procedure. In addition to 
those waves, we identify in our analysis new gravity waves 
throughout the analyzed domain. In particular, note the 
high energy density (dark red) regions near the 50 km, 
300 km and 400 km altitudes all with m* w 0.25. These 
constitute new waves. 

Second, we identify a region consistent with a turbopause 
~50 km thick, centered at approximately 400 km altitude. 
This is the region where the width of the sub-spectrum con- 
taining significant energy increases markedly. In this re- 
gion, molecular diffusion becomes comparable to eddy diffu- 
sion and gravity waves — growing in amplitude as they prop- 
agate upward — break, transferring energy into the higher 
wavenumbers. Above this region, energy is lower across the 
entire spectrum and the spectrum itself is much steeper (i.e., 
narrower) than that for the turbopause region. Above the 
turbopause region, the atmosphere becomes inhomogeneous, 
separating out into layers of different molecular species. Past 
studies have placed Jupiter's turbopause at the ~ 5/ibar 
[Festou et al, 1981] and ~0. 5/ibar [Yelle at al., 1996] levels, 
based on observations and modeling. Our result supports 
the latter location. 

Third, in the lower part of the analyzed domain there ap- 
pears to be a ducting region, a region with a sharp jump in 
N. As already noted, such a region can serve as a source 
of gravity waves. Horizontally propagating gravity waves in 
this region have previously been observed (see, e.g., Fig. 1), 
which have been suggested as waves trapped in a "leaky" 
duct. The ducted wave travels horizontally by undergoing 
internal reflections at the boundaries; however, part of the 
wave escapes the duct to propagate vertically. Because the 
wavenumber with maximum energy in the £ map can be 
traced down to the ducting region, it is likely that the waves 
have come from there. 

Note that above the ducting region the number of local 
peaks in the spectrum generally reduces with altitude (the 



centroid of the spectrum at each height is shifted to lower 
wavenumber). Also, the magnitude of the peak energy in 
the low wavenumber (white line in Fig. 3) is high at first, 
then decreases, and then increases again along the white 
line, until the topside of the turbopause region at ~ 425 km 
above the 1 bar level. This is indicative of wave satura- 
tion or encounters with a critical layer, where 1 = locally, 
for high number wavenumbers as they propagate upward. 
The second or third multi-peaked altitude regions, between 
~ 200 km and ~ 300 km above the 1 bar level, may be due to 
wave breaking and secondary wave generation from breaking 
layers. 

As alluded to earlier, m* differs from the required 
wavenumber m. The latter is the wavenumber that would 
be observed by a probe travelling in a vertical direction. 
However, throughout most of its entry phase, the probe had 
a shallow angle of attack (~7° below the horizontal), which 
changes the wavenumber observed by the probe during its 
passage in the upper part of the analyzed domain. We use 
the geometry of the probe's path to obtain the true vertical 
wavenumber, 

/ tan 7 tan /3 \ » 

m = - m , (3) 

y 1 - tan 7 tan B J w 

where 7 is the angle of attack and /3 is the angle the wavevec- 
tor makes with the horizontal. We make no correction for 
the relative motion of the probe with respect to the wave 
since the probe's velocity is supersonic (indeed, hypersonic, 
with up to Mach 51) for much of the entry phase. Gravity 
wave phase speeds are subsonic. 

Although 7 = 'y(z) is known from the probe's trajectory 
[Seiff et al, 1998], /3 is not. To estimate this quantity, we 
assume that the probe's trajectory is vertical in the period 
just before parachute deployment. This is not far from the 
actual situation since the probe's angle of attack was 83° just 
before the parachute was deployed, near the lower boundary 
of our domain. In this region, m* « m; and, since 

8 = arccos | — ) , (4) 

\Vk 2 + m 2 J V ' 

we have B « 85°. Also, in this region m « 0.25 km -1 (giv- 
ing vertical wavelength of w 25 km) , in good agreement with 
previously recovered values [Magalhaes et al, 2002; Arregi 
et al, 2009]. Now, all the parameters required to recover 
/ have been obtained, and a vertical profile for / can be 
estimated. 

Finally, to obtain uo from I a value for the horizontal 
phase speed c is required. This is not well known. How- 
ever, equatorial gravity waves with a phase speed of about 
100 m s~ greater than the winds have been observed in im- 
ages of Jupiter's clouds returned by New Horizons [Reuter 
et al, 2007]. This gives c « 180 m s _1 . The Doppler Wind 
Experiment (DWE) on the probe has measured the zonal 
winds to be approximately this speed in the deeper atmo- 
sphere [Atkinson et al., 1998], suggesting that the waves 
could have been driven by convective overshoot below the 
clouds. This value also agrees well with the 140-195 m s _1 
range reported for the similar region, at the probe entry 
site [Magalhaes et al., 2002]. Somewhat smaller values (ap- 
proximately 70-145 m s _1 ) have been obtained for waves at 
different longitudes [Arregi et al, 2009]. Such variation in 
c does not change the shape of the uo profile we recover, 
but the magnitude would be reduced. Variations of c with 
altitude would increase the uncertainty as well, but they are 
expected to be small and not fundamentally change the pro- 
file. This gives the profile for uo{z) shown in Fig. 4a. Here 
we have used c = 180 m s^ 1 . The standard deviation shown 
indicates the variation in the profile given by the various 
smoothing window sizes considered in constructing £ . 
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4. Discussion 

The DWE reported flow speeds that increase with depth 
reaching 170 m s _1 at the 5 bar level and then remaining 
high at lower levels [Atkinson et al., 1998]. The upper part 
of the DWE profile is shown in Fig. 4a. The wind speed at 
the bottom of our profile agrees well with that at the top of 
the DWE profile. Our profile shows increasing zonal wind 
speed with altitude (up to ~100 km above the 1 bar level). 
This implies that the wind speed is a minimum near the 
cloud-top level. This is similar to what has been observed 
from studies using the thermal wind equation [Flasar et al, 
2004; Simon-Miller at al., 2006]. However, in our profile this 
high speed wind is not a jet, as the zonal wind speed does 
not diminish appreciably with altitude near this level. In 
fact, the speed is roughly constant throughout the strato- 
sphere, starting from this level. There are some fluctuations 
of the order of 20 m s _1 in the wind speed in the thermo- 
sphere. The Richardson number, Ri = N 2 (du /dz)~ 2 , is References 
greater than 1/4 for the entire profile, which indicates the 
flow is stable with respect to Kelvin-Helmholtz instability. 

Temporal variations in the temperature profile of 
Jupiter's equatorial stratosphere have been observed to have 
a period of 4 to 5 years and are therefore known as the quasi- 
quadrennial oscillation (QQO) [Leovy et al, 1991; Friedson, 
1999]. These oscillations have been linked to variations in 
the zonal wind observed via cloud-tracking on Jupiter with a 
period of ~4.4 years [Simon-Miller at al, 2007]. Further, a 
jet located just north of the equator, derived using the ther- 
mal wind equation applied to data gathered in December 
2000 and January 2001, has been linked to the QQO [Flasar 
et al., 2004]. A jet of similar magnitude, just south of the 
equator, has been observed in observations gathered in 1979 
[Simon-Miller at al, 2006] (which fits well with a period 
of 4.4 years), suggesting this jet may be oscillating about 
the equator. No such jet is visible in our profile. However, 
the probe entered Jupiter's atmosphere 5 years before the 
northern jet was observed. Thus, our profile is consistent 
with such an oscillation as the jet would be located south of 
the equator, away from the probe entry site, at the time the 
probe entered. 

Gravity wave creation and dissipation is associated with 
heating and cooling of the atmosphere. Ignoring molecu- 
lar viscosity, which is only important above the turbopause, 
the heating rate is proportional to the vertical gradient of 
energy flux — specifically, 



explore the troposphere of Jupiter to depths of 100 bar or 
more. This will provide better insight to the source and 
behaviour of gravity waves in the troposphere, possibly al- 
lowing better limits to be derived for the wave-vector angle /3 
in our analysis and better understanding of the mechanisms 
that generate the gravity waves. The planned Jupiter Icy 
Moon Explorer (JUICE) mission will directly study gravity 
wave activity and zonal winds in the stratosphere of Jupiter, 
extending the result present here. 

Acknowledgments. C.W. is supported by the Science and 
Technology Facilities Council (STFC), and J.Y-K.C. is supported 
by the STFC PP/E001858/1 grant. The authors acknowledge 
H. Thrastarson, I. Polichtchouk, and C. Agnor for useful discus- 
sions. We thank the reviewers for helpful suggestions. 



~dt 



1 dF z 

pc p dz 



(5) 



We estimate the energy flux, F z = p4>w, from A9 using 
the polarization relations to derive the perturbations in the 
geopotential <j> and the vertical velocity to. The flux comes 
from the zonal average of the product of these quantities. 
Since the product is not available to us we average over a 
single vertical wavelength instead. The heating rate profile 
thus obtained is shown in Fig 4b. It can be seen that heat- 
ing is small through the stratosphere. Indeed, on average 
it is very close to zero. It is only in the thermosphere that 
the average heating rate deviates from zero with a magni- 
tude of ~ 0.5 K per Jupiter rotation. The flux responsible 
for heating in this region is small (F z < 10 -3 W m~ 2 ), 
compared to radiation fluxes measured in the troposphere 
(3 < F z < 16 W m~ 2 ) [Sromovsky et al, 1998]. Never- 
theless, the peak heating rate magnitude is about 50 times 
larger, partly because the lower thermsopshere is much less 
dense than the troposphere. 

The circulation of Jupiter's stratosphere is important for 
understanding the planet's circulation as a whole. The lo- 
cation of the turbopause is essential for understanding the 
coupling between Jupiter's upper atmosphere and the cir- 
culation in the lower atmosphere. The Juno mission will 
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Figure 1. Gravity wave in Jupiter's atmosphere (Voy- 
ager image 16316.34). This wave is propagating horizon- 
tally in the troposphere. Many such waves are captured 
in the Voyager images [Flasar and Gierasch, 1986]. Im- 
age contrast has been enhanced to improve the visibility 
of the wave. The wave's location is indicated by the ar- 
row. Image courtesy NASA/JPL-Caltech. 
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Figure 2. Relative potential temperture perturba- 
tion and Brunt-Vaisala (buoyancy) frequency profiles 
for Jupiter's atmosphere. All altitudes are relative to 
the 1 bar pressure level. The profiles are based on 
data gathered from t = —173.055 s to t = 218.170 s; 
t = s is the point of parachute deployment. Our anal- 
ysis used the complete set of acceleration data collected 
from both accelerometers with the exception of one out- 
lier at t = —157.742 s. The acceleration data collected 
by the probe includes a small spurious oscillationSeijff 
et al. [1998] which we did not smooth out as its effects 
on the calculation was negligible in the atmospheric re- 
gion we analyzed. The bottom 25 km of the profile is 
from direct measurement of the temperature with the 
first 15 s of direct temperature measurements removed as 
these were anomalously high, a, The vertical profile of 
potential temperature perturbations, AO, scaled by the 
background value, 8. The perturbations show wavelike 
oscillations throughout the stratosphere and lower ther- 
mosphere. Note that the short wavelength oscillations 
in the layer between 25 km and 50 km are due to accel- 
erations caused by buffeting of the probe as its velocity 
became subsonic, b, The Brunt-Vaisala, frequency is the 
maximum frequency that a gravity wave can have and 
shows the stability of the atmosphere against convective 
stability. 
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Figure 3. Moving Lomb-Scargle periodogram of the po- 
tential temperature perturbation from Fig. 2a. Regions 
of high spectral energy are in red while low energy is in 
blue. The variation of the observed vertical wavenum- 
ber m* is shown (white line). The line is produced by 
identifying local energy maxima and constructing a line 
joining them, avoiding local minima, starting at the level 
of the duct in the troposphere. Note that there are other 
waves within the periodogram we do not consider in our 
analysis. One such wave, the region of high energy at 
m* w 0.61 km - between 180 km and 210 km altitude 
(labelled Wl) has been previously identified as a saturat- 
ing gravity wave Young et al. [2005]. Another region (la- 
belled W2) has also been previously identified Young et al. 
[2005]. The region at an altitude of around 400 km shows 
a broadening of the range of wavenumbers with increased 
spectral energy. This is indicative of the turbopause, the 
region where waves break and the atmosphere begins to 
become heterogeneous. 
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Figure 4. The vertical profiles of the background zonal 
wind speed Mo and heating rate, dT/dt). a, The zonal 
wind is shown (solid) with the average standard devia- 
tion of the variation across smoothing and periodogram 
windows indicated. The zonal speed profile found by 
the Doppler Wind Experiment is shown at the bottom 
(dashed) for comparison, b, The vertical profile of the 
heating rate is shown (dotted), as temperature change 
per Jovian rotation (9.925 h), with the average stan- 
dard deviation of the variation across smoothing and pe- 
riodogram windows indicated. 



